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Abstract 

The conductivity of a finite temperature 1+1 dimensional fermion gas described 
by the massive Thirring model is shown to be related to the retarded propagator of 
the dual boson sine-Gordon model. Duality provides a natural resummation which 
resolves infra-red problems, and the boson propagator can be related to the fermion 
gas at non-zero temperature and chemical potential or density. In addition, at high 
temperatures, we can apply a dimensional reduction technique to find resummed 
closed expressions for the boson self-energy and relate them to the fermion conduc- 
tivity. Particular attention is paid to the discussion of analytic continuation. The 
resummation implicit in duality provides a powerful alternative to the standard 
diagrammatic evaluation of transport coefficients at finite temperature. 



1 Introduction and motivation 

The evaluation of transport coefficients at high temperatures in terms of Feynman dia- 
grams in a weakly coupled theory is a subtle and highly involved task [1, 2]. Firstly, they 
are often proportional to the mean free path of the scattering processes, which increases 
as coupling strength decreases. Further, in the relevant limits of vanishing external mo- 
menta and energy, higher loop diagrams can be as important as lower loop diagrams if 
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they are sufficiently infrared sensitive. The failure of perturbation theory that both of 
these observations imply requires careful and clever resummation of diagrams. 

In lower dimensions we know that this resummation is a reflection of the fact that 
the relevant degrees of freedom are not those of the quasi-particles, but of dual degrees 
of freedom. This is well understood in the case of Luttinger liquids, electrons close to the 
Fermi surface in one dimension (quantum wires). Such a system [3, 4, 5] has a dual rep- 
resentation of its charge modes in terms of free bosonic fields, which provide the relevant 
degrees of freedom. This leads to a huge simplification in the calculation of conductivity 
[6], not easily visible (if at all) from electron Feynman diagrams. The predictions that 
this simple duality permits have been confirmed in experimental measurements [7] of the 
conductance of GaAs quantum wires, although the finiteness of the system enforces a 
modification [8] of the naive picture. 

In this work we explore the advantages of using this fermion-boson duality when 
calculating conductivity for relativistic quantum fields in 1 + 1 dimensions. The hope is 
that the infinite-order non-perturbative resummation of the quasi-particle modes can be 
replaced by a few terms in the series for the dual degrees of freedom i.e. that duality does 
the resummation for us. Further, the usual calculation of transport coefficients relies on 
linear response theory, and the simplifications implicit in the dual resummation suggest 
that we can go beyond linear response, otherwise impossible. 

As our example, we will concentrate on perhaps the simplest non-trivial theory dis- 
playing conductivity: the Massive Thirring (MT) fermion theory, with Minkowski La- 
grangian density 

Cmt$, V>] m )V - l^M^x), (1.1) 

where = (p,j) = is the fermion number current. 

This is dual to the sine-Gordon (SG) boson model, which is described by 



£s G [0] = ^0^-^cos(A0). (1.2) 



provided the renormalised coupling constants are identified as 

A 2 1 



4tt 1 + g 2 /n 



1.3) 



a 
A 2 



pm, (1.4) 



where p is the renormalisation scale and m the renormalised fermion mass. 

We note that, with multiplicative renormalisation, if m = 0, then the massless 
Thirring model is dual to the free bosonic field. If we introduce a chemical potential for 
the fermion field then, in the limit of no antiparticles, we recover the Luttinger limit. 
However, in QFT the massless limit is unnatural for quarks or electrons and we are 
obliged to be more sophisticated. 

While 1+1 dimensions may seem unrealistic in relativistic quantum field theory, the 
links of our model with conformal field theory means that these results are relevant to 
more exotic situations. For instance the study of decay of quantum normal modes of 
a classical field in the presence of a black hole, all in AdS space, can be linked through 
AdS / CFT correspondence to linear response in a conformal field theory in 1+1 dimensions 
[9]. 

Further, as a precursor to understanding the deconfinement transition in QCD in 3+1 
dimensions, the two-dimensional fermionic Thirring model can be put in correspondence 
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with a compact U(l) gauge theory in four dimensions [10], by virtue of Parisi-Sourlas 
dimensional reduction [11]. From this viewpoint the chiral bilinears of the MT theory 
correspond to the monopole-antimonopole pairs of the four-dimensional theory, and the 
monopole condensation that signals confinement has its counterpart in the chirally broken 
phase of the MT theory [12], that we discussed in detail in an earlier exploration of this 
duality [13]. 

In [14, 15] it has been shown that the SG/MT models are also equivalent at finite 
temperature T and fermion chemical potential /i. Here, we will consider the fermion 
response to an electric external field (conductivity) as a working example of the use 
of duality and dimensional reduction in the calculation of transport coefficients. This 
amounts to an evaluation of the full retarded boson propagator. In this way, the self- 
energy of the bosons is directly related to the conductivity and charge screening of the 
fermions. Furthermore, we will show how to relate the dynamical information contained 
in the boson self-energy, calculated in the imaginary-time (IT) formalism, to the static 
properties of the fermion gas. In addition, we will use dimensional reduction techniques 
to analyse the boson propagator. Dimensional reduction is meant to be valid at high 
temperatures compared to the fermion mass, strong fermion couplings and large distances 
compared to the inverse temperature. This regime has been shown to be very useful in 
[13], where it was used to resum exactly the pressure and the fermion condensate at finite 
T and ji. 

The structure of the paper is the following. In section 2 we derive the relationship 
between the fermion induced currents in linear response theory and the SG retarded 
propagator. The retarded propagator can be read off from the IT one by analytic con- 
tinuation. We devote section 3 to analyse the imaginary-time SG propagator. First, we 
write it as an infinite sum of mass insertions and then we explore its relevant proper- 
ties and the relationship with the fermion gas at finite chemical potential. In section 
4 we discuss the high temperature limit and the dimensional reduction regime for the 
SG propagator. Analytic continuation from imaginary energy to real energy is a fraught 
exercise, and this model is simple enough to show how subtle one needs to be. This 
will become clear in sections 4 and 5 where we will relate the analytic continuation of 
the high T propagator with the physical conductivity. We will show that one needs to 
impose physical conditions on the propagator in order to have a physically meaningful 
answer for the transport coefficients. Many details of the calculation have been collected 
in the Appendices. Appendix A contains some useful results about thermal propagators 
and analytic continuation, while Appendix B contains some details of the calculations 
performed for the SG propagator. 

2 Induced fermion currents and the boson retarded 
propagator 

In practice, the only transport coefficient that we will study is electrical conductivity, a, 
the response of electric currents to an electric field, since it is for this that duality can be 
easily exploited. In our case, the electric field will be an external classical field appearing 
as a c-number source in our action. It is convenient to absorb the gauge coupling e into 
the gauge fields so that the covariant derivatives are = <9 M — iA^ and the gauge field 
has dimension one. The interaction term added to the MT Lagrangian density is then of 
the form 

5V(t)=Ux,t)A^(x,t), (2.1) 
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where the external classical electro-magnetic potential is and j M = (p,j) = is 
the fermion number current. Note that using these conventions, the gauge field kinetic 
term is of the form — F^F^ /(4e 2 ) so the electric charge, e, has dimensions of energy. As 
a gauge choice, we will take Af = in what follows, so that the external electric field is 

E d {x,t) = -(d/dx)A$(x,t). (2.2) 

with dimension two. The electric field satisfies Maxwell's equations in 1+1 dimensions, 
which read 

e 2 j c i(x,t) = -E d (x,t) 

2 / ^ dE cl (x,t) 
e p d (x,t) = — . (2.3) 

With these conventions the conductivity a is defined by 

j = aE, (2.4) 

where j is the current density and E is the electric field. The plasma is isotropic and 
homogeneous so the conductivity is a simple scalar. 

In defining conductivity, we are implying that the response of the plasma to an electric 
field is simply linear in that field. This weak field limit is calculated in quantum field 
theory using the theory of linear response (see [16, 17]). In terms of the additional 
interaction term, 5V, of (2.1), the terms linear in electric field come from the first term 
of a perturbative expansion in powers of SV. From linear response theory this is then: 

roc roc 

({6Wj ll (x,t))) = -i / dt' dx'A^x',t')(([j,(x,t),j (x',t')]))9(t-t'). (2.5) 

J to J-oo 

The Sj^x, t) is the difference between the currents with and without the external pertur- 
bation. The 6^ indicates that it is the term linear in SV, and hence linear in the electric 
field. ((■)) stands for a thermal average taken with respect to an equilibrium density 
matrix based on the unperturbed Hamiltonian given at time earlier than to G K. The 
perturbation, SV is assumed to be switched on only after t and we are looking at the 
induced current perturbation at time t G K > to- 

For the MT model of (1.1), any thermal average involving the fermion current j M 
can be replaced by a boson thermal average in the sine-Gordon model of (1.2). Duality 
means the fermion current has an exact representation in terms of the scalar SG field 
alone, namely 

jp{x,t) = ^e^0(x,t), (2.6) 

with A defined through (1.3). The proof of the above statement can be found in [15]. If 
the fermion fields interact electromagnetically, then we would identify g = 0(e) in the 
effective theory described by (1.1). We feel no need to make this identification, and treat 
g as an independent parameter to be varied irrespective of e. 
Putting this together for the current then gives 

(\ \ 2 poo POD o o 

— J jf dt' J dx'A^x^tyit-t'^^mx^Aix^t')]}}. 

(2.7) 

Using integration by parts on the x' variable and using (2.2) gives an expression in terms 
of the electric field. A resulting boundary term at spatial infinity can be ignored in 
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practical problems. The time derivative can be made to act on both expectation value 
and theta function by using the equal-time commutation relation [(f>(x,t),(f>(x',t)] = 0. 
This then leaves us with the linear response of the fermionic current in terms of the 
applied electric field given in terms of the full real-time retarded propagator A R (x,t) of 
the sine- Gordon scalar field 



(\ \ 2 poo poo o 

A J jf dt'j dx'E d (x',t')^A R (x-x',t-t'), (2.8) 

where 

A R (x-x',t-t') = e(t-t f )(([ ( f ) (x,t), ( j ) (x',t')]}} = | + °°^ f + °° ' ^e-**e*"A R (u,,p). 

(2.9) 

Taking to to minus infinity (ignoring for now possible subtleties in this limit [18, 19]) 
allows the Fourier transform to be taken, giving 

((S {1) 3i(w,p))) = (J^j E d (u;,p)ujA R (uj,p). (2.10) 

The above result is nothing but Kubo's formula [17] for the SG/MT system. Thus, the 
conductivity as defined in (2.4) is given in terms of the SG retarded propagator as (note 
that ((j))=0 without external fields) 

A x 2 



<r(w,p) = luA r (lu,p). (2.11) 

Proceeding exactly in the same way for the zero component of the current yields the 
induced charge density: 

((5 {1) jo(w,p))) = E d {u,p)pA R {u,p). (2.12) 

So we see that all we need for the linear response to an electric field is the full retarded 
sine-Gordon scalar propagator at real Minkowski energies, which we shall analyse in detail 
in next sections. 
Note that: 

• The induced current is conserved, d^S^j^ = as it should. 

• In the free case, the retarded propagator is T-independent (see Appendix A) and 
so are the fermion conductivity and induced charge density. 

We conclude this section by considering briefly the case when the bosonic degrees of 
freedom are those of a free field. 



2.1 Free bosonic modes 

We see from the duality conditions that the massless (mo = 0) Thirring model is dual to 
the massless and free bosonic field, for all couplings g 2 . In this case the conductivity is 
given by 

a(u,p)=(^j coA^cp), (2.13) 
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where A° R (u,p) is the retarded propagator for the massless free field. 

In this, and subsequent sections, it is convenient to keep a small nonzero mass /i 
for the scalar field and in the end we will take — > + . We shall discuss retarded 
propagators in great detail later. For the moment it is sufficient to quote the result 

A° R (u,p) = — — - 2- ( 2 - 14 ) 

{00 + it) 2 -v -i4 

We recover the conductivity from (2.13) as 

A \ loo 



a(uj,p)=[ — \ f , (2.15) 

\2-k J (LU + Id) 2 - p 2 - /ig 

In particular, the real part of the conductivity for constant applied field is (after taking 
/i ^0+) 

A\ 2 _ 1 



(7(0;, 0) = tt ^- j 5(u) = IT77 - 5(u). (2.16) 

We shall see in section 5 that, if we apply a static field only to a finite part of the system, 
then the conductance is 

G= K^) = 27r(l+\v 7 r)- (2 ' 17) 

Suppose we had not appreciated that the massless Thirring model was dual to a free 
bosonic theory. We would then have attempted to calculate the conductivity (or conduc- 
tance) as a series expansion in g, using the bilinear fermionic forms for the directly 
in (2.5). The first term in the interaction picture expansion is the simple one-loop term, 
which gives [20] 

a(u, 0) = 6(u) or G = —. (2.18) 

2ir 

We now see the power of duality in resumming the series in the fermion coupling 
constant g. Further, the reduction in the conductance due to the presence of repulsive 
interactions (g 2 > 0) has an exact counterpart in Luttinger liquids [8]. 

In subsequent sections we see how duality aids resummation for the massive fermion 
theory, through the dual sine-Gordon theory. 



3 The sine-Gordon propagator in imaginary time 

We will adopt a similar approach to that of [15, 13] in expanding in fermion mass about 
our results for the massless fermionic theory and its dual bosonic free theory counterpart. 

As before, it is necessary that the bosonic calculations are moderated by a small 
nonzero mass /xq- In [15, 13] it is shown that observables such as the pressure and the 
quark condensate are /i -independent and hence infrared finite. This is true also for 
correlators evaluated at different space-time points [15, 13] as long as they involve at 
least one derivative of the scalar field. The conductivity as it stands in (2.11) should be 
also Ho independent, because it is proportional to the time derivative of the two-point 
function, and so is pAr(uj,p) in (2.12). 

The imaginary time propagator is obtained by differentiating twice the generating 
functional 



A t (x,t) = ——— -^Z SG \J-T\ 

Z SG (T) bJ{x,T) bJ(0) 

6 



j=o 



Z SG [J;T] = Nfs f d<j> exp{- / dr f [C SG [<t>] + J(x)<f>(x)} \ , (3.1) 

J periodic L Jo J —oo ) 

where r G [0,/3], /3 = T _1 and Zsg{T) = Zsg[0;T] is the SG partition function, which 
coincides with the MT model one Zsg(T) = Zmt(T) as showed in [14, 15]. 

The SG generating functional (3.1) can be given explicitly as a power series in a (see 
equation (3.16) in [15]). With a/A 2 = mp, this is an expansion in fermion mass about 
the free boson theory. Setting J = one has the SG partition function: 



Z SG (J = 0,T) 1 + 

71=1 

2n j-1 



Z B (T) 

["dr, f^dx) (UUmxj-x^rj-r^r^) , (3.2) 

yj=l J ° J -°° J \j=2k=l J 

where Zq(T) is the free boson partition function (equal to Zq(T), the free fermion 
partition function in 1+1 dimensions), and 




€j -^ - j j = n + l, n ..,2n ■ (3 ' 3) 



The Q 2 function is given by 



Q 2 (x, r) = sinh(-^ '-) sinh(-^ '-), (3.4) 

and A°(i,t) is the free boson IT propagator in the /i — ► + limit [14, 15]: 

A° (x, t) = ~ In pl(3 2 Q 2 (x,T). (3.5) 

It can be useful to visualise expressions like (3.2) in terms of diagrams representing 
the terms in a expansions such as (3.2) or (3.12). Each expression has 2n integrations 
over 2n space-time coordinates, each represented by a vertex. They split into two types 
according to the choice (3.3). The n vertices associated with €j = +1 and integrals over 
the (xi,Ti) to (x n ,T n ) coordinates we will denote by a closed circle. The remaining n 
vertices have ej = —1 and integrals over the (x n+ i, r n+1 ) to (x 2n , T 2n ) coordinates, and we 
use an open circle. It is a crucial rule that there are always as many open as closed circle 
vertices as a result of a superselection rule, the derivation of which was a key element of 
[15]. This rule is nothing but taking the no — ■> + limit so that physical quantities are 
IR finite. There is also a factor of a/A 2 = mp per vertex. Note that in terms of the 
fundamental field, tag corrections (A(0) factors) of the field vertices (when thinking 
in terms of an expansion of the cos(A0)) have been absorbed leading to a (T/p) A / 47r 
renormalisation factor for each a, e.g. through the use of thermal normal ordering [15]. 
Thus, each of these vertices represents a factor of 

I dT > L "'if- a(T) := a b) = aa h?) • (3 - 6) 

with the e factors associated with each vertex modifying the "propagators" attached to 
them. The second form for a(T) is in terms of the bare coupling and the UV cutoff A, 
as given in [15]. 
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In every diagram each of the 2n vertices are connected once to every other vertex by a 
"propagator", drawn as a double line. The double line linking the j-th and k-th vertices 
represents a factor of (Q x / 27r ) e<e J so that lines connecting one open and one closed vertices 
contribute to the denominator while those connecting vertices of the same type are part 
of the numerator. Without the factors of e from the vertices, we have 

[Q 2 (x,r)]- AV4 * = ^=~exp[A 2 A„] 

~ + \x<Q + lA«Q + ... (3.7) 

In terms of the field each Q propagator line represents many different diagrams 
with vertices of all even orders coming from the expansion of the cos factors in the 
Lagrangian. However a A expansion is not a good way to appreciate the result as a 
whole because of IR problems. In fact, note that the original expression for the partition 
function (and the same will happen for the propagator) was fio independent, while every 
single diagram in (3.7) depends on fiQ, giving rise to divergent terms as fio — ► + , which 
is the limit where the identification Q 2 ~ exp[— 47rA ] done in (3.7) is valid. We will 
therefore keep the double line resumed propagator and we will discuss in section 4 its 
interpretation in the high temperature limit. 

Since the Q functions (3.4) are invariant under Xj 4=> x^ then no sense of direction 
need be assigned to these double line propagators, and there is a symmetry between 
diagrams related by an exchange of open and closed vertices. 

Summing up, the vacuum diagrams can be represented by sums over diagrams of the 
type illustrated 




Note that the interaction was based on a cos(A0) and not a cos(A0) — 1 factor, so there 
is a physically irrelevant shift in the action and resulting partition function. Also note 
that though this is a calculation of the partition function Z and not In Z, the diagrams are 
connected. Finally, we remark that the above diagrams are in position, not momentum, 
space and there are only 2n — 1 independent integration variables in each term in the 
sum due to translation invariance, as it is emphasized in Appendix B for the case of the 
propagator. 
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Following similar steps as in [15] for the partition function, one can write the full IT 
propagator as an infinite series in the renormalised a (3.6). Again the /io — > + IR limit 
must be dealt with properly, but in the end we can express the result in the same terms 
as the partition function with a couple of extra rules. First, as in normal perturbation 
theory, one will 'pull out' two free propagators ~ ln(— Q 2 /4ti) (see Appendix B for 
details) of (3.5), which we will denote with a single line. These connect to any of the open 
and closed vertices appearing in the partition function, but for each additional single 
field line attached to a vertex j one multiplies by a factor of —iejX/ (47r). However, a key 
result in this model is that one has only one or two free propagators contributing to the 
full propagator. One does not get chains of self-energy insertions separated by a single 
free propagator. We will comment about this observation in section 4. 

Thus, the full propagator has only four types of contribution, a single free propagator 
and three types of diagrams with interactions. 



A T (x,r) = A°(*,r) + |^£rfW), 



n=l 



rf n \x,r) = r^(x,r)+r^(x,r)+r^(x,r). 



,(2n) 



(2n) 



,(2n) 



(3.10) 
(3.11) 



The three types of interacting terms are given in terms of a similar formula, differing only 
in how the two free propagators are convolved with the Q 2 propagators. Writing those 
terms as {AA} A (A = ±, ±=p, ±±) the general expression is of the form (see details of 
the derivation in Appendix B): 



if;W) 



A 

47T 



a (T 



A 2 /4vr' 



-\ 2n 



2n j-1 

x ( aa )a i n n - ^ t j - 

\j=2 k=l 




dxj 



(3.12) 



where 



{AA} ± 

{AA} ±T 
{AA} ±± 



2nA^(x 1 , t 1 )A^(x 1 - x, n - r), 
-2n 2 A° (xi, ri)A° (x 2n ~ x, r 2n - r), 
2n(n - l)A° T (x!, t^A ^ -x,t 2 -t). 



(3.13) 
(3.14) 
(3.15) 



Note that we have exploited the space-time translation invariance of an equilibrium sys- 
tem as well as the symmetry between open and closed vertices to choose specific space- 
time variables for the free propagators with the remaining equivalent choices being ac- 
counted for by simple combinatorial prefactors. Such factors are obvious when we write 
these in our diagrammatic notation as 



(3.16) 




+ ^ — , +(•< — ► ") + ... (3.17) 




T^l(x,r) = / \ +(• — o) (3.18) 

+ (. < — > o) + ... (3.19) 




rf±l(x,T) = / \ +(.< — ► o) » X 02 *^ +(. < — >°) + ..(3.20) 



The last two diagrams above can be combined in interesting ways, exploiting the fact 
that they differ only by the exchange of an open for a closed vertex on an external leg, 
with a compensating change on an internal vertex. 

It is not difficult to check that all the integrals appearing in the above expressions 
for the propagator are finite in the infra-red, i.e, for large spatial distances, following 
the same arguments as in [15] for the partition function. In addition (see Appendix B 
for details) the only /^-dependence in (3.12) appears in the free part A°(x, r), in the 
/jl — > + limit. In Appendix B it is also shown that the integrals in (3.12) are UV (short 
distances) finite provided A 2 < An. We shall restrict here to A 2 < An, where the theory 
is superrenormalisable (see comments in [15, 13] in this respect). On the other hand, the 
free propagator is UV divergent. Therefore, all the IR and UV divergences are in the 
free part of the full IT propagator for A 2 < An. Hence, every term in the a-expansion in 
(3.12) is infra-red (and UV) finite, which indicates that this is the appropriate expansion 
we should look at to avoid the problems mentioned in the introduction, and related here 
to the naive A expansion. 

We also remark that (3.12) is independent of the scale p, since the explicit dependence 
is exactly compensated by the implicit scale dependence in a (see [15] for more details). 
A conventional choice of scale is p = m, the mass of the fermion. This is particularly 
interesting at high temperatures T ^> m where the a expansion in (3.12), originally 
seen as an expansion in fermion mass m 2 , becomes effectively a high temperature series 
in (m 2 /T 2 ) 1-A / 8?r [13] (see section 4). Thus, at high temperatures, we can consistently 
truncate the series (3.17), (3.19) and (3.20), although we shall not do so yet. Notice that 
we fix the mass of the fermion, meaning in particular that a — > + in the A — > + limit. 
In other words, even in that limit we do not recover a massive free boson theory, as one 
would naively expect from the lagrangian (1.2), by expanding to lowest order in A. 

Now let us Fourier transform the propagator (3.12) to momentum space A(iu n ,p) 
with discrete frequencies cu n . We can write the result as: 

A T {iu n ,p) = A^(iu n ,p) [l + i: T (iuj n ,p)A^(iuj n ,p)] , (3.21) 
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where A° (icu n ,p) is the free propagator in momentum space: 



A T (iu n) p) 



1 



< + p + n 



2 ' 



(3.22) 



where we have kept the dependence with the small fi® for reasons to become clear below. 

Note that S is not a self-energy in the usual sense, but rather it is just the truncated 
full propagator. It is 1PI but the full propagator is not an infinite sum of such insertions. 
The conventional Schwinger-Dyson form for the propagator in momentum space reads 
schematically A T = A^(l — II T A) where IT T is the true self-energy. Therefore, the 
imaginary-time self-energy is related to S T as 



U T (iuj n ,p) 



T, T (iuj n ,p) 



1 + ^ T (iu n , p)A° T (iu n) p) 



(3.23) 



In more detail, the non-trivial part, the T, T (iuj n ,p), can be written as (see details of 
the calculation in Appendix B) : 



E T (iuj n ,p) 



A 2 Z»(T) 
^Z SG (T)^\n\ 



^2 ( n \ 



a 

2>7 



A 2 /4tt' 



2n 



E^(iuj n ,p), 



(3.24) 



with 



2" p(i r+oo 
j=2 ^° l<k<j<2n 



l=k+l 



l=k+l 



€j€ k \ 2 /An 



(3.25) 



Note that, from (3.25) it is clear that £^(0,0) = and Sy(iw m ,p) = iuj m ,p) = 
T^TiioJm, —p). Besides, from the asymptotic behaviour of the Q 2 variables as x'j — > ±oo 
it is not difficult to see that the Taylor series of Y*T{iu m ,p) around p = for fixed oj m is 
well defined, so that we can write 



^ T {to rn = Q,p) = Y j a k (T)p 2k . 



(3.26) 



k=i 



In addition, note that, from (B.5) one has [S T A^](0,0) = (since X)fc=i e k = ^)- 
To make this result compatible with (3.26), one has to keep the small mass in the free 
propagator /i ^ till the very end of the calculation, since 



[E T A°]K t = 0,p) 



a lP 2 + 0{p 4 

p 2 + 14 



(3.27) 



so that [E T A^](tu m = 0, 0) = if fi ^ 0, and so we will do in the following. 

Here we will be mainly interested in the analytic continuation of E^u^p) to real 
frequencies u and its behaviour for small frequencies and long wavelengths (small p). 
Before discussing the analytic continuation, we will derive an interesting relationship 
between the lowest order coefficient a\ in the momentum expansion (3.26) and the MT 
model at finite density. It is a good example of how the use of duality can yield interesting 
and unexpected connections between the boson self-energy and the fermion gas. 
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3.1 The lowest order boson self-energy and the fermion charge 
density 

From (3.26) and (3.25) we have 



«?(T) = ^n/ dr> n 

x <( n(n- l)[x' 2 ] 2 -n 2 



l<fc<j<2n 


" 2n 








J=2 . 





Q 2 ( E ^ E ^) 



=fc+i 



=fc+i 



(3.28) 



Here, a?(T) denotes the C?(p 2 ) coefficient of SJ(0,p). Let us now perform back the 
change of variables (B.7)-(B.8) and put the system on a finite length L, so that dxj — > 

I~L/2 dx r ^ n ^ ne en d we snan take -L — > oo. Thus, (3.28) becomes 



2n 



a l( T ) = -Tr\[ / dT 3 I dx 3 II [Q 2 ( X 3 ~Xk,Tj -n)] 

x [n(n - 1) (x 2 - xif - n 2 (x 2n - %i) 2 } ■ 



£je k \ 2 /<lir 



(3.29) 



At this point, let us recall that the MT model partition function at nonzero chemical 
potential /i can be written as [15, 13] 



Z MT (T,fj) = Z F (T)exp 
with 



PL 



2(tT + < ? 2 ) 



oo / \ 2 



n=0 



a (T 



2A 2 VP 



A 2 /4tt' 



2n 



F 2n (T,/i),(3.30) 



2 ™ Z"' 3 f n/Z A 2 /4 

F 2n {T,fi) = Y[ / drj / tZajj JJ [<2 2 (^ - x k ,Tj - T k )] e3tk ^ exp 



»L/2 

" L / 2 l<fc<j<2n 



Zq(T) = Zq(T) being the free fermion (or boson) partition function. 
Consider now 



^2 z ™ 
E e J x J 

(3.31) 



:^2n(T,/i) 



/A 2 \ 2 2n ^ f L/2 
= ~ 2 \az) 11/ dT 'J I dx i II [Q 2 (^i -^> r i -^)] 

^=0 V 47r / j=1 -/0 J-L/2 i< fe<j < 2 „ 

x {ra; 2 + n(n — l)xix 2 — n 2 xix 2n } , 



ej€ fe A 2 /47r 



(3.32) 



where we have relabelled the Xj variables as explained in Appendix B for similar calcu- 
lations. Further relabelling allows us to replace in the above integral: 



nx\ + n(n — l)x±x 2 — n 2 xix 2r , 



1 n 2 

-n(l - n) (x\ + xj- 2x x x 2 ) + — (x\ + x\ n - 2rr 1 x 2n ) 



= ~\ [ n ( n - l)( x 2 - xi) 2 - n 2 (x 2n - xi) 2 ] , 
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(3.33) 



and therefore, by comparing with (3.29) we find 



47T / Att 



PL\X 2 J dp? 



2 d 2 



F2n{T, fX) 



(3.34) 



M=0 



Let us turn this relationship into one involving physical observables. The pressure of 
the MT model gas is 

P M t(T, fi) = \im —In Z MT (T,fi), (3.35) 



and the fermion charge density is 



Pmt (T, p) = —P MT (T, p) = -^- 2 + p c (T, p), 
op 7r + g 2 



(3.36) 



where the first term in the r.h.s. is the massless Thirring model charge density [21] and 
the second one is given in [13]. Then, collecting our previous results and recalling that 
Z MT (T,p = 0) = Z SG (T) [15] we find for the lowest order boson self-energy: 



16tt 2 d 



A 2 dp 



(3.37) 



This equation provides an exact relationship between the leading low-energy boson self- 
energy and the fermion charge density subtracting the massless part, as a direct conse- 
quence of the fermion-boson thermal duality. 



4 High temperature approximate propagator 
4.1 Static observables 

Before looking at dynamical quantities at high temperature, such as conductivities, it 
is worthwhile recalling how one can study time-independent quantities in the SG/MT 
system. Just as in other models in other space-time dimensions, the large temperature 
means that the Euclidean time dimension becomes very small and one expects dimen- 
sional reduction (DR) to occur [22, 23, 24]. 

For the SG/MT model, the nature of dimensional reduction and its application to 
static quantities was studied in [13]. The key observation was that the Q 2 variables of 
(3.4), the essential building blocks of all the exact SG/MT expressions, have a large- 
distance limit of 

Q 2 (x,T)^\e W 2 -^ for V M < 1 « M. (4.1) 

If the relevant distance scales for the physics of interest are much longer than /3, then we 
will be able to use this approximation for the Q 2 . We will refer to this limit as DR, since 
the dependence of Q 2 on imaginary time r disappears. For the static quantities, such 
as the pressure, the precise conditions required for DR to be a good approximation were 
shown in [13] to be that either T ^> m with g 2 > (A 2 < Att) or T > m with g 2 ^> 7r 
(A Att), with m held fixed. Here, m is the fermion mass at the scale p = m. 

The main utility of DR for the SG/MT system is that it allows us to write a-expansions 
such as (3.30)-(3.31) in terms of a classical one-dimensional gas of 2n charged particles 
whose positions on the line are labeled by the Xi, subject to the Coulomb potential 
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V(xi,x 2 ) oc \xi — x 2 \. This system was solved exactly a long time ago [25, 26], so that 
in the DR regime we can resum the a series and obtain exact results for the thermal 
SG/MT system [13]. In particular the pressure is a static quantity and thus the phases of 
SG/MT can be obtained. As shown in [13], one can identify a "molecular" phase, where 
fermion condensates tend to pair forming "molecules" (in analogy with the behaviour of 
Coulomb charges in [25, 26]) which are responsible for the chiral symmetry restoration 
as T — > oo. There is also a lower T regime, or "plasma" phase, where condensates pair 
less easily and the chiral symmetry is broken. 

For quantities with no space-time dependence, such as the pressure, the applicability 
of DR is intuitively obvious and can be made mathematically precise as in [13]. On 
the other hand, as it is also shown in [13], if we are dealing with space-dependent but 
static objects, such as correlators of fields separated in space but not time, then the DR 
approximation for Q 2 is a good one provided we work in the limits T ^> m, A <C 4n 
and we enforce an extra condition. That is that any external spatial distance scale, x, 
satisfies \x\ ^> /3/ir. In terms of Fourier components this means that the approximations 
can only describe 'soft' physics, that is physics of spatial momenta p < T. 

Of course, in the DR regime the simple form for the Q 2 (4.1) means we can obtain 
closed results for the different expressions. However here we wish to study time dependent 
quantities, and it is not at all clear that a time- independent approximation for the Q 2 , 
such as the DR form, can have anything to say about the study of dynamical quantities. 
In fact, as we will see now, the time- independent approach amounts to neglect all but 
the n = Matsubara mode. However, the n ^ modes do play a crucial role in physical 
situations. 

As first pointed out by Pisarski [27] the study of dynamical quantities at high temper- 
ature in four-dimensional weak coupled theories requires that one resums an infinite set of 
diagrams as the high temperature can compensate for the small coupling. This has been 
since been developed in detail from several perspectives [28]. The original diagrammatic 
viewpoint of Braaten and Pisarski [29] is that one has to resum HTL (Hard Thermal 
Loops) — the leading temperature dependent contributions from diagrams of arbitrary 
order in coupling. This leading temperature dependence comes only from the hard modes, 
where energy and momentum is at least of order T. However we can also view this as 
producing an effective action describing the 'soft' physics in 3+1 dimensions, that is on 
energy or momentum scales of order E,p < gT or less, where g is a gauge coupling or 
equivalent. To get such an action, the hard modes must be integrated out of the theory. 
The point is that such an effective theory can describe dynamics of soft processes not 
just static quantities. Moreover, the HTL effective action has non-trivial dependence on 
energy (at least for relativistic fermion and gauge fields in 3+1 dimensions) and so differs 
from the energy independent DR effective action. Though soft-momentum dependence 
is described by the HTL effective action, it is only for static quantities that it reduces to 
the DR action. The results of [28], if not the general ideas, are specific to 3+1 dimensions 
and to relativistic theories. One of our tasks here is to study this question in our models. 
We will show in next sections that one has a hierarchy of soft and hard scales similar to 
the HTL one and also emerging from an effective resummation. 

4.2 The imaginary- time SG propagator at high T 

We will work in the following limits: T > m (m is the fixed fermion mass), p <C 2nT 
and A 2 <C 47r for the reasons explained above. 

In the T m limit, it is justified to keep only the 0(a 2 ) term in the a-expansion for 
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the IT S T in (3.24). We follow similar arguments as in [13]. That is, we shift r- — > fir'j 
and x'j — >■ (3x'j in (3.25) so that we can write (3.24) for p = m as 



Z T {iu; n ,p) = T 



2 A 2 Z B (T) ~ /J. 



^Z SG (T) 



n=l 



1 / 2\ 1-A 2 /8tt' 

1 / wr \ ' 

2 1 T2 J 



2n 



f%.(2nin,p/T), (4.2) 



where is a dimensionless function obtained from in (3.12) by taking f3 — 1, so 
that it depends weakly on T for small p/T (it is independent of T for p = 0). If we 
do the same with the partition function Zsg(T) as discussed in [13], (4.2) yields a well- 
defined expansion in m/T (remember that A 2 /4?t < 1). The advantage of using this 
expansion instead of the naive A expansion is that it is free of the problems mentioned in 
the introduction, like the bad IR behaviour. 

Hence, taking the leading order in the high-T expansion gives: 



j: T (iuj n ,p) 



A 2 /m 2 \ 1_A2/4w 



[h(iu n ,p)-h(0,0)], 



with 



h{iuj 



p(3 poo 

iu n ,p)= dr dx e WnT e- tpx [Q\x,t)\ 

JO J-oo 



-A 2 /4tt 



(4.3) 



(4.4) 



Taking into account that the full IT propagator is A T = A T + S T [A T ] 2 , keeping only 
the term (4.3) can be interpreted diagrammatically as: 



+ 2 



(4.5) 



with all other diagrams ignored. 

The next step is to analyse the high-T limit of h(iu> n ,p). For that purpose, we will 
consider the DR for the Q's discussed in [13]. That is, we have: 

Q 2 (x, t) = ]- [cosh(27rTx) - cos(2ttTt)] , (4.6) 

so that 

[Q 2 (x,T)]~ X2/4n = 2^-^1*1/2 [l +e -4 7 rWT_ 2e -2 7 rNT cos(27rTr) ]-A 2 /4 7 r_ ^ 

Note that for large distances \x\ — > oo, the integrand of (4.4) is screened with the 
thermal mass scale 

m T = A 2 T/2, (4.8) 

which ensures that the function h is IR finite and in the end will allow us to obtain IR 
meaningful results for the conductivity. 
If we expand in (4.7) 



[1 + 4 



-A 2 /4tt 



= E 



r [i - a 2 /(4^)] 



fc=0 



r [1 - k - A 2 /(4tt)] T(k + 1) 



A 



(4.9) 
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where A = e" 4 ^l T - 2e- 2 ^l T cos(2ttTt), we get 



h(iw n ,p) = XUZS 5, 




(4.10) 



r [1 - k - A 2 /(4tt)] p 2 + (2nT(2k - /) + m T ) 2 ' 



Our approach will consist in keeping only the dominant terms in the DR limit (p <C 
2nT and A 2 <C An) for the sums in (4.10) but keeping the exponential term exp(— 7tit\x\) 
in (4.7), in order to reproduce the correct IR behaviour. For instance, consider the n = 
contribution to (4.10). It is not difficult to see that the dominant contribution is the term 
k = I = m = in the above sum, namely, 

h(0,p)^2 X2 ^^^. (4.11) 
p z + 

while the contributions with k > 1 are subdominant. For instance, for k = 1 we have the 
non-leading contribution 



which is subdominant compared to (4.11) in the limits p <C 2nT, A 2 <C An. In turn, we 
see why it is important to take also the small A limit (in the sense explained above) as we 
had anticipated in the previous section. That is, has to be considered as a quantity 
of the same order as p and both m T ,p <C 2nT. In the language of HTL, p and m T are 
soft scales and T is a hard scale. Thus, we see why it is consistent to keep m T without 
expanding further in A in the exponential in (4.7). 

We can do the integral (4.4) numerically and check our previous approximation. For 
n = 0, before doing the x-integral in (4.4) we have in the DR regime 



The comparison between the l.h.s of (4.13) numerically integrated and the asymptotic 
expression on the r.h.s is showed in Figure 1. The approximation is worse as X 2 /An 
increases or Tx decreases, as expected. Nevertheless, for values of A even as large as 
A = 3 it works remarkably well for |x| > \/T. In practice it is important that the 
condition A 2 An does not have to be enforced rigidly, since small A requires large g 2 . 
Although (2.16) is valid for all g 2 , we would like our result to be more than a super-strong 
coupling result at high temperature. 

The above procedure can be followed also for the n^O modes, for which the leading 
order term in the sum is k — I — \n\ (note that n^O implies k ^ and hence I < k < 2k). 
Expanding also to leading order in A 2 /47r we get 




(4.12) 




(4.13) 




(4.14) 



(4.15) 
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Figure 1: The dashed line is the l.h.s of (4.13) evaluated numerically and the solid line is 
the asymptotic expansion in the r.h.s, for two different values of A. 

where lo^ = 2nkT. 

An important remark is in order here. The expressions (4.14) and (4.15) are equivalent 
in the DR limit. That is, they differ in 0(p 2 /(2nT) 2 ), 0(m^/(2nT) 2 ), <D(\ 2 /4it) terms. 
However, that does not mean that their analytic continuations to real energies will be 
equally close for all energies. For that reason we have given them separately because they 
provide a clear example of how the analytic continuation of slightly varying imaginary- 
time results (including arbitrary non-leading order terms) can give different physical 
answers. As a matter of fact, and for reasons to become clear below, we will also introduce 
a further modification of (4.15): 

»""i^k * § s <r > s - © 

(4.16) 

where m s = ©(my) = 0(X 2 T) is a soft scale mass. We could take m s = m T , although 
(4.16) allows for more general situations (see below). 

Finally, combining our results in this section with those in section 3.1, we get a high-T 
expression for the slope of the MT charge density at the origin simply by replacing in 
(3.37) the coefficient a± of the p 2 term in £y(0,p), which we readily obtain from the 
zero mode contribution (4.11). The expression thus obtained is consistent with the result 
found in [13] for the charge density. 

4.3 Analytic continuation to real energies 

We now have an expression for the SG imaginary-time propagator Ay(iu; n ,p) based on 
h(iuj n ,p) of (4.14)-(4.16). However, all the physical information we will need for dynamical 
quantities such as the conductivity is encoded in the retarded SG propagator at real 
energies. So one needs to analytically continue from the values Ax(iuj n ,p) at discrete 
imaginary frequencies iu n to arbitrary complex uj. 

The work of Baym and Mermin [30] showed, as explained in appendix A, that by 
enforcing a particular set of conditions, those given in (A. 20), one produces a generalised 
propagator, A(u n ,p), a function of complex energies u, whose values at u = E ±ie with 
E real, give the retarded and advanced propagators, e.g, eq.(A.15). 
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Let us try to apply the procedure of Baym and Mermin to our expression (3.21) for 
the IT propagator, and try to find the retarded equivalent of (3.21). The generalised 
propagator will have the form 

A{u,p) = A°(u;,p) + iX(u;,p)[A (u;,p)] 2 . (4.17) 

The analytic continuation of the first term in (3.21), a free IT propagator, is well known 
and used as an example in Appendix A and in section 2.1. The generalised continued 
function £(o>, p) can either be defined from (4.17) as the second term in the full propagator 
divided (truncated) by two generalised free propagators, A . With a little work one can 
see that £ is indeed the analytic continuation of the IT function £y, satisfying a slightly 
altered set of BM conditions (A. 20) as its large energy behaviour can be such that £c<j~ 4 
vanishes for large \u\. It is then a short step from this generalised function to the retarded 
version using (A. 15), and we will find the form 

A R (E,p) = A R (E,p)+iZ R (E,p) [A° R (E,p)] 2 , (4.18) 
where the first term is simply the free retarded propagator 

and the second term we can write as a retarded function multiplied by two free retarded 
propagators. Physically, we will be interested only in the small p behaviour of T,(u,p) at 
high temperatures. 



4.3.1 Approximate analytic continuations 

In the previous sections, we have analysed the IT propagator to leading order in the high 
T expansion (4.2). In fact, to that order, and according to (3.23), £ r rs —U t , with 
£^ m (4-3) and the imaginary-time self-energy. One must bear in mind though that 
the analytic continuation of perturbative approximations to (3.23) for arbitrary complex 
energies must be done carefully, taking into account the poles of the free propagator and 
the £ functions. We will discuss this issue in section 4.3.3 below. 

We have been able to find approximate expressions for h(icu n ,p) in the DR limit in 
(4.14)-(4.16). However, as discussed in greater generality in Appendix A, it is possible 
to find different functions whose analytic continuations all give the same IT result ap- 
proximately, i.e. up to non-leading order terms in the DR limit. In this section, we will 
give explicit examples of such functions for the case of h(icu n ,p). All of them satisfy the 
conditions explained in Appendix A, i.e, such functions h(cu,p) are analytic off the real 
axis, huj~ A vanish for large |a>| and h(u> = iuj n ,p) coincides with (4.14)-(4.16) to leading 
order in DR. In section 4.3.2, we will show that additional conditions may be imposed to 
ensure that the analytic continuation gives unique and meaningful physical answers. Let 
us then consider the following cases: 

1. The first choice we could make is to neglect everything but the zero mode contri- 
bution in (4. 14)- (4. 16), which is dominant in DR with respect to the n^O modes, 
i.e, we take just h(iu n ,p) = h(0,p)5 n o. As discussed in Appendix A (see (A. 24) and 
(A. 25)), analytic continuation gives a function h\{uj,p) which vanishes everywhere 
off the real axis. Therefore, for E EM., 



h 1R (E,p) = h 1 (E + ie) = 0, 
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(4.20) 



so that only the constant term h(0, 0) would survive in the analytically continued 
retarded propagator of (4.3). In this regard, we note that when only the leading 

zero mode contribution is taken, i.e, when one replaces just [Q 2 {x, r)] A ^ w by 
2A /27r e -A t\x\/2 j n (4.7) ; one can obtain closed expressions for the full self-energy 
sum in (3.24)-(3.25) up to 0(p 2 ) for any u n , in terms of the Coulomb gas pressure 
and charge density, following similar steps as in [13]. However, as the example of 
the free propagator in Appendix A shows, this analytically continued propagator 
may be a very bad approximation (we will be more precise below). More realistic 
dynamics requires knowledge of the contribution of the n^O modes as well, even if 
each heavy mode is of a smaller order in our approximations than those neglected 
in calculating our zero mode contribution. 

2. With this in mind, for the second case we take (4.14) as the starting point. From 
this form, one can apply the BM conditions (A. 20) and derive a unique analytic con- 
tinuation to real energies. In this case, as so often, a straight iu n — > u replacement 
in the functional form is essentially sufficient and we find 



h 2 (u,p) = - 2 " /2nX [ s : = e [i m{uj)] _ e [-ImH] (4.21) 

p A + [m T — %suj\ 

Note that the because of the modulus on the integer k in (4.14), the form of the 
generalised function changes in the upper and lower half planes. This is common and 
through (A. 15), this then corresponds to the retarded and advanced functions being 
distinct. Thus, the retarded function of (4.21) for real E is h,2R,(E,p) = h 2 (E + it) 
with s = 1. We remark that the result (4.21) obeys the BM conditions (A. 20), 
reproducing all the Matsubara modes in (4.14) including n = (see the discussion 
in Appendix A). 

3. If we take (4.15) instead of (4.14) as our IT function, then we would have, as our 
analytic continuation, 

_2A 2 /2ti-^2 

which again gives hs(u,p) ~ h{iuj ni p) up to non-leading order terms in DR. 

4. Finally, consider the fourth expression, obtained by analytically continuing (4.16): 

K{u,p) = -X 2 {^) (P+ 2 m %^ } 2 ■ (4.23) 



Clearly, we could obtain infinite variety in our analytic continuations just by modifying 
the original high-T imaginary time expression by non-leading terms. The important point, 
as emphasized in Appendix A, is that the difference in physical quantities obtained by 
analytically continuing functions differing at non-leading order should be also of non- 
leading order. The four examples we have shown above are enough to show that this is not 
necessarily the case, unless further restrictions are applied, as we will see below. In fact, 
note that the analytic structure of the retarded version of those functions is very different. 
For instance, while h 2 in (4.21) for s — 1 has complex single poles at uj — ±p — irriT, 
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Figure 2: The solid line is h^(u,p) as given by (4.22) and the dots represent the IT values 
at u = iu n calculated numerically from the definition (4.4), for n < 5. 

/i3 in (4.22) has real poles at u — iE p with E p = \Jp 2 + m 2 -, and in (4.23) has 
branch cuts for real u 2 > p 2 + m 2 in addition to the singular behaviour at oj 2 = E 2 . In 
general, it is not difficult to realise that by moving the pole of the AC function in the 
complex plane by soft amounts, but preserving the zero mode contribution, the IT values 
at n ^ remain unchanged to leading order. This can be easily achieved by replacing for 
instance [uj 2 — p 2 — m 2 ^] by [uj 2 — p 2 — m% + uM + (u,p)} for Imu > in the denominator 
of h, with M + <C 2ttT a regular complex soft function such as the denominator does not 
vanish for Imu > and so on for Imu < with another function M_(cj). Thus, while the 
u = iu n values are approximately equal, the contributions near the respective poles are 
arbitrarily different. The poles of the retarded propagator are related to physical modes 
at T ^ 0, or quasi-particles, while the branch cuts of the self-energy have to do with their 
decay rate. However, as it is clear from the previous discussion, we cannot determine the 
presence of soft poles unambiguously. There is no novelty in this as a general problem. 
For example, the whole programme of Pade approximants is based on how best to exploit 
such ambiguities. We will address this question again in section 4.3.3 although in the 
end the physical results will not be affected by our ignorance about the precise analytic 
structure. 

Before proceeding, let us comment that one can check numerically that the above 
functions indeed match the Matsubara modes approximately. For instance, in Figure 2 
we have plotted the discrete n-values against the function h 3 (u,p) in (4.22). We see that 
the agreement is quite good, given the numerical uncertainties, for small p and small A, as 
expected. It gets worse as A or p increase, although it still gives a reasonable agreement 
for p < 2tcT. Note that as p increases, the n ^ mode become of the same importance as 
the zero mode, but the magnitude of the latter becomes much smaller. Recall that if we 
plotted the approximation (4.20), it would be just a function that vanishes everywhere 
except at the origin, where it matches the zero mode. From Figure 2 we see that this 
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is a cruder approximation than h 3 in the range u 2 < (2irT) 2 . This gives support to the 
idea, explained in Appendix A, that we need to combine the information about the n^O 
modes, essential for AC, with that on the one soft mode n — 0, even if we can only do 
this approximately. 



4.3.2 A physical condition 

We are interested in the fermion conductivity, which as we have seen is related to the 
SG retarded propagator. In particular, as we will see in section 5 (see eq.(5.2)) the time 
evolution of the conductivity is directly related to the Fourier transform of the retarded 
propagator: 

/°° dE 
^e- iEt A R (E,p), (4.24) 

A R being given in (4.18), where in the high T limit, 

A 2 /m 2 \ 1_A2/47r ~ 
^p)^-m 2 T 2 [^) h R (E,p), (4.25) 

with 

h R (E,p) = h R (E,p)-h(0,0), 

where h(0, 0) ~ 2 a2 ^ 2n A 2 / 'm\ to leading order. 

Consider now the contribution of the second term, i.e the term, to (4.24). As 
discussed in Appendix A, the retarded function h R (z,p) can be also analytically continued 
to complex z, simply starting from the retarded h R (E,p) for real E and replacing E by 
z. This function, by construction, has its poles and branch cuts in the lower half plane. 
On the other hand, the generalized analytic functions h we have been discussing in the 
previous sections satisfy Schwartz's reflection principle h(z*) = h*(z) and the property 
h(z) = h(—z). This means that h R (—z,p) = h* R (z*,p). 

Therefore, taking t > and closing the integration contour from below the real axis, 
we find the following contribution to A R (t,p) given by the poles at uo 2 = oo 2 = p 2 + /i 2 ,: 



1 d f -2vri 



27T<9/iQ [ 2uj p 



h R (u p ,p)e luJpt - h R (-u p ,p)e 



d 1 



iuipt 



1 



Reh R (p,p) (sin pi — ptcospt) — lmh R (p,p) (cospt — ptsinpt)j 



(4.26) 



where /i'^ = (d/du)h R (uj,p). Note that /ir does not depend on /z (see our comments 
about the IR behaviour in section 3) and that the expression (4.26) is real, as it should, 
since it contributes to the conductivity. Therefore, from (4.26) we see that A R (t,p) has a 
term growing linearly in time (and so does the conductivity) unless h R (p,p) = (although 
h' R (p,p) will be in general different from zero). 
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Hence, if we want the conductivity to remain bounded in time, as seems physically 
reasonable 1 this gives us the following condition for the retarded £ function on the mass 
shell of the SG massless field: 



Z R (E 2 =p 2 ,p) = (4.27) 

for e — > + . This is an extra, physical, condition, in addition to the BM mathematical 
ones discussed in Appendix A. 

From the IR behaviour of h R as p — > + , discussed in previous sections, it is not 
difficult to see that once we demand (4.27), the rest of the contributions in (4.26) remain 
bounded in time when taking the Fourier transform also in the p variable. As for the 
contributions of the poles of h R in the lower half plane, say at uj = ±E(p) — ij(p) with 
7(p) > 0, their contribution is damped at long time as exp(— jt) if 7 7^ (as it is the 
case for instance for the retarded version of (4.21) where 7 = m T ) and one can check 
that for 7 = 0, although there is no exponential damping, the contributions to A R (t,p) 
remain bounded as well, since one gets sin Et and cos Et instead of smou p t and cosu p t as 
in (4.26), so that all the /xq dependence is in the t-independent contributions. 

Finally, if h R (z,p) has also a branch cut, as it is the case for h iR (z,p) = h 4 (z + ie,p) 
with hi in (4.23) for z = —ie ± E' with E' > \Jp 2 + m 2 , the branch cut contribution to 
A R (t,p) gives also a bounded contribution, as long as h R (E,p) is integrable for E > E cut 
where E cut is the cut endpoint. For instance, for the case of h± R with m s = m T , such a 
contribution is proportional to 



~iE- 1 Shl < £t > 



Ie, ""(B 2 -p 2 ) 2 (£;2_ B |)i-g' 

with Ep = \/p 2 + rriip, which is perfectly finite and bounded in time, since the integrand 
behaves like (E — E p )~ r with < r < 1 near E p . The same is true for m s 7^ m T , as it is 
discussed in section 5, see eq. (5.11). 

Of the four cases considered in the previous section, only /i 3 and h 4 satisfy the on-shell 
condition (4.27) and they yield respectively: 



^sr(E,p) ~ -fir 



E 4R (E,p) w -ji T 



{E + tef-p 2 
(E + ie) 2 — p 2 — m 



2 ' 

T 



m 2 T m; x2/27T [p 2 + m 2 -(E + xeff 1 ^ 
(E + ie) 2 —p 2 — m\ 



+ 1 



with 



on 1-A 2 /4tt 



A = 2^ 2 /-m 2 . (4.28) 

Note that p^ m 2 in DR, but we have not made any assumption about the relation 
between m and tut and therefore between \i T and m T . However, it is natural to take m 
and m T of the same order, since both are soft scales much smaller than T, and therefore 
p\ m\ . In fact, most of our results below will be expressed as corrections of order 
/iy/m^. This means that the temperatures should be at least as high as T>m/ X 2 but they 
can be even much higher, which would have simply the effect of making the perturbative 
corrections much smaller. 



1 In addition, an unbounded correction would become eventually larger than the leading order, so that 
the high-T expansion would be meaningless. 
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4.3.3 Analytic structure and quasi-particles 

Here, we will discuss some relevant features of the retarded propagators we have just 
obtained. The poles of the retarded propagator in the lower half plane give the dispersion 
law for quasi-particles in the thermal bath [17], which are the dual SG bosons in this 
context. In moving from the MT model, with two fundamental degrees of freedom, to 
the sine-Gordon model with one bosonic degree of freedom, we appear to have lost a 
degree of freedom. However, the second bosonic degree of freedom was trivial and was 
integrated out when exploiting the duality [15]. For instance, for m/T — > 0, bosonic and 
fermionic modes both contribute the same amount to the free energy in 1+1 dimensions. 
We must take account of this trivial bosonic mode in the thermodynamic quantities, but 
for dynamic quantities, only the single SG mode is necessary for equivalence with the two 
massive Thirring modes. 

On the other hand, the branch cuts of the self-energy U(E,p) along the real axis 
correspond to the decay rates of those quasi-particles and the corresponding discontinuity 
U(E + it) — U(E — ie) = 2TmfI across the cuts can in principle be calculated by cutting 
the self-energy diagrams. Cutting rules and discontinuities in the self-energy at finite 
temperature have been extensively analysed in the literature [31, 32, 33]. To one loop 
in perturbation theory [31] and for the case of a single particle with mass m, the self- 
energy has branch cuts at s = Am 2 , 9m 2 , ... (s = E 2 — p 2 ) corresponding to the decay 
into two-particle, three-particle states and so on. These are the T = cuts although 
the discontinuity across the cuts is T-dependent. New "thermal" cuts may appear, even 
to one loop, corresponding to processes of emission and absorption of particles in the 
thermal bath, not allowed at T = 0. 

In our case, even to leading order in the high T limit, we are dealing with a non- 
perturbative resumation of diagrams in the coupling constant A, as for instance in (4.4). 
As commented above, this is essential in order to obtain a meaningful IR behaviour for 
physical quantities such as the transport coefficients and it has forced us to use analytic 
continuations to approximate imaginary-time results. Thus, we do not have a clear inter- 
pretation in terms of Feynman diagrams that can be cut and therefore we do not know a 
priori what should be the analytic structure of the self-energy on the real axis. To make 
matters worse, and as emphasized in previous sections, with our DR approximation we 
cannot determine the position of the self-energy branch points or the poles of the retarded 
propagator. One can move the poles by soft amounts (i.e, of 0{m T )) and the analytic 
continuation will still match the IT values to leading order in the DR high-T limit. 

The point we want to make here is that, despite this apparent limitation, we can give 
meaningful physical predictions. The key point is that, once the extra condition (4.27) is 
imposed, the errors produced by changing the poles or the branch cuts are perturbatively 
small. We will check this explicitly in section 5. Before that, let us discuss in some more 
detail the analytic structure of the AC propagators. First, let us define the retarded 
self-energy in terms of the retarded propagator as customary, 

A*(w, p) = \ r. (4.29) 

(Lu + ie) 2 -p 2 -Tl R {uj,p) 

Note that we cannot write (4.18) as A R ~ i(u 2 — p 2 + S_r) _1 since, as explained in 
section 3, we do not get an infinite series of S insertions. Remember also that the poles 
of A R (uj,p) lie in the lower half of the complex plane uj and the retarded propagator for 
real energies E is A R (E,p). 

Consider first (4.21), which gives for (j, — 
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.//f, [s + 2i(u; -He)m T ] + s [s - mf, + 2i(u + ie)m T ] , . OAN 

A 2j r(u;,pJ = « Tl 2 i t/ 1 ' (4.JUJ 

s 2 [s — + 2z(u + ?e)mTj 

where s = (a; + ie) 2 — p 2 and /if, defined in (4.28). 

We see that A 2 r has a double pole at s = and single poles at c<j = ±p — im T . On the 
other hand, the retarded self-energy defined through (4.29) satisfies lmUn(E + ie,p) ^ 
for any E real different from zero, so that it has a cut along the whole real axis. 
Note also that Ur(uj,p) is singular for the complex oj solution of /if, [s + 2i(u + ie)mr\ + 
s [s — mf, + 2i(u; + ie)mT\ = 0. 

We compare this with the solution (4.22), which gives 

A 3R (u,,p) = i S + t %-^ . (4.31) 
s[s — 

Now, A 3 r has single (and real) poles at s = and s = mf, and the self-energy is real 
on the real axis except at the singular point {oj + ie) 2 = p 2 + mf, — /if,. 

The analytic structure of the above two solutions is completely different. As we have 
discussed already, the position of the poles and branch cuts change within DR. A different 
story though is the physical interpretation of the above results. First, note that the two 
solutions above for the retarded propagator have a common feature: the self-energy is 
singular close to the new "thermal" poles (i.e, those different from s = 0). That is, the 
distance between the singular points and the "thermal" poles in the oj plane is propor- 
tional to ht, which is much smaller than the positions of the poles themselves, which is 
0{jtit)- The values for which the self-energy diverges are related to pinching or end point 
singularities [33] and it is therefore not clear that those poles should be interpreted as 
quasi-particles energies or damping rates. However, the behaviour of the two above solu- 
tions near the pole at s = is very different. Thus, while A 3R (s — > 0) ~ i{\ — /if, / mf,) 
is a massless pole with a T-dependent residue, A 2 r(s — > 0) ~ i(±2ip/if,) / '(±2ip — m T )s~ 2 
is a double pole. A double pole is an indication of the breakdown of the Schwinger-Dyson 
expansion around the massless propagator [33] and in fact we have seen in the previous 
section that it leads to unphysical behaviour. Note that the condition (4.27), which is 
satisfied by A 3R but not by A 2 r, ensures that the pole at s = is a single pole. There- 
fore, our physically acceptable solution for the retarded propagator is consistent with a 
massless dispersion law. However, we want to make clear that our approach does not 
allow to say whether the quasi-particles dispersion law really remains massless or there 
is a "thermal" true pole of 0(rriT) as it happens in the HTL approach. All that we 
can say is that m^ plays the role of a screening mass in the sense explained in previous 
sections. In addition, note that the contribution to physical observables from the pole of 
A 3R at s = mf, is reduced by a (/if,/mf,) factor with respect to that of the s = pole, 
so that the thermal contributions in DR are always perturbatively small with respect 
to the leading order. Another indication of the ambiguities in the interpretation of the 
retarded propagator is the result (4.23), which also satisfies the condition (4.27) and we 
could have chosen to have a branch cut starting at m s = 2m T , as expected from a true 
thermal quasi-particle pole of mass m T [33]. 

Summarizing, the DR approach does not allow us to fix uniquely the position of poles 
and branch cuts of the retarded propagator. Although the results suggest the appearance 
of a thermal mass scale mr, it is not clear within our approximation scheme whether 
this is the mass of the quasi-particles. Thus we see that the double-lined "propagator" 
introduced in section 3 really has a pole but this does not simply correspond to a stable 
particle of mass m T . The important point is that the physics we are interested in, i.e, 



24 



the conductivity, does not depend on our choice of functions to be continued analytically, 
once the physical condition (4.27) is imposed. For instance, taking (4.22) or (4.23) will 
produce the same answer, to the order we are considering here. 

5 Physical applications. The conductivity. 

Transport coefficients are usually discussed in Fourier space and the expressions such as 
(2.11) for conductivity, used to motivate the calculations of the full scalar propagator, 
were in energy and momentum. Our interest is in the DC conductivity, i.e. the u,p — > 
limit of (2.11). However, the free field propagators Ao appearing in expressions such as 
(2.11) are poorly defined in this limit because they are massless. The DC conductivity 
then seems to depend on how the u,p — > limit is approached although, in simple cases, 
as in (2.16), there may not be any difficulty in making a choice. For more complicated 
expressions more care is needed. Such unsatisfactory behaviour is well known in thermal 
field theory for zero-momentum Green functions [18, 19], often coming from Landau 
damping cuts across zero energy rather than a massless pole as here. 

Our solution is to remain in coordinate space where one can specify a realistic ex- 
perimental situation, with systems of finite size examined for finite times. The lack of 
analyticity found at zero momentum and energy at finite temperature can be associated 
with space-time causality. After all, Sj should be zero if we are looking at a point that 
is space-like separated from a non-zero electric field. The usual formulae will then give 
a — 0. On the other hand for the same problem at a time-like separated point, a will 
be non-zero. Likewise a simple implementation of a DC conductivity measurement, the 
E,p — > limit, implies that a constant electric field has been applied for all times. If 
the system has finite correlation times, this may be an acceptable approximation to an 
experiment where the fields were set up a 'long time' before. However, the presence of 
poles or cuts, at zero energy and momentum, suggests that there are long time scales in 
the problem, so one can not simply assume that such a DC current can be switched on 
adiabatically in the distant past. In practice, working in time and space is not especially 
difficult, and it is much more physical, so it greatly simplifies the physical interpretation 
of real-time calculations in finite temperature problems. 

Turning to our model, we must first specify the electric field. From Gauss' law the 
electric field due to a single charge in vacuo is constant in size, merely switching sign at 
the location of the charge, i.e. E oc q(9(x — Xi) — 9(x 1 — x)) is the field due to single 
static charge at x±. A constant field over a finite region x\ to x 2 comes from having two 
charges of opposite signs at X\ and x 2 . Thus we choose 

E d (t, x) = E9(t - t ) [9(x - x^ - 6(x - x 2 )\ , (5.1) 

where E is a constant and x 2 > X\ for simplicity, so that E = E if x\ < x < x 2 , and is 
zero otherwise. This is then a close analogue of the usual large parallel plate capacitor 
problem encountered in 3+1 electrostatics. 

However, we are turning on the field suddenly, so the two charges appear instanta- 
neously at time to- The field is switched on over a region which is not initially time-like 
separated, so we must be careful later with the interpretation. This means that we have 
to calculate 

\ 2 POO P%2 p ,-72 

6jW(X) = i^jf dt'J dx> J ^e-^^')e^').(-ip )A il (po,p) 
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_ A " p V „ f d2 P r -ipo(t-t ) ip(s-»,Q A fl(P0»P) /. 2 x 

" Ax** 2-f*> 7 (27r)2 6 6 p ' ^ 

where si = +1, S2 = —1 and X = (x,t). 

As shown in section 4, the full retarded propagator is made up of two parts, as given 
in (4.18). The first part is the free contribution and for the second one we have obtained 
explicit expressions in the DR high-T limit. Let us analyse those two contributions 
separately. 

5.1 Free boson term 

Consider the first term, the single free propagator (4.19) and its contribution 5ji 1 \x) to 
the total 5j {1) (X). This gives, for t > t : 

(1) A 2 -^ f°° dp sin [p(x - Xj)} sm.{up{t - t )) 

5ja (x) = ^ h 3 L * — v ^ — 

J"=l,2 



where u p = \Jp 2 + /j 2 , and in the second line we have taken the /xo — > + limit after 
performing the p integral. This result is clearly IR and UV finite, even if /x — > + . 
However, it is also causal, i.e. the current is zero if its space-like separated from the 
electric field region. The regions where the functions are non-zero are shown in Figure 
3. The 8{ is one in region I - the region inside the forward light cones from the point 
(to,Xj), i.e. the space-time point where the charge at position j was switched on. It is 
zero elsewhere. The Q\ is then one only in region II to the right (positive x) of this and 9 J 3 
is one in region III to the left (negative x). Each of the terms in the sum for j represent 
the effects of switching on at to the electric field of a single charge at Xj. Note that 
the regions II and III are space-like connected to (t ,Xj) yet the change in the current 
is non-zero in these regions. This is because we have actually switched on the electric 
field of these charges instantaneously, and the static field is non-zero everywhere. A more 
realistic experiment would be to switch on the charges slowly (charging up a parallel plate 
capacitor). However this behaviour is not important for the case at hand. In fact, when 
we add the two terms in the sum for j together we get 



where 



5£\X) = j atmax [6 n + 2(t-t )L~ 1 6 23 + [(t-t )-(x-x 2 )]L~ 1 6 21 

+ [(x-x 1 ) + (t-t o )]L- 1 13 ], (5.4) 
A 2 

ja,max = ~^2~EL, (5-5) 



e kl = e\el l = x 2 - Xi . (5.6) 



The regions where the 6m functions are non-zero are shown in Figure 4. If the distance 
between the charges is L, i.e. the electric field E is on over a region of length L, then the 
current profile rises until at time t — t — L/2 > it reaches a state where there is a central 
region where the current saturates at j a ,max and is proportional to L for constant electric 
field strength E. This region extend a length t — t — L/2 either side of the mid-point 
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of the electric field region for t > L. The current profile then drops linearly to zero at a 
distance t — t + L/2 either side of the mid-point. Its then zero beyond this, as it must 
as these regions are not in causal contact with the region of electric field (i.e, the regions 
22 and 33 in Figure 4). This is shown in Figure 5. 

Thus there are two conclusions. After an initial rise time, the current in the region 
of electric field is constant proportional to L. Thus the contribution to the conductivity 
from the first term in the full propagator is 



a = 



8tt 2 ' 



(5.7) 



as in (2.16). Thus the current is always proportional to the voltage, EL, and independent 
of L so the conductance (the inverse of the resistance) is constant, at G — A 2 /87r 2 , as 
in (2.17), whatever the length of the material is studied! However, life is not quite so 
simple. Note that in fact the current is also nonzero outside the region of the electric field. 
In fact, the current is reaching the same constant level everywhere as fast as causality 
allows. This is nothing but the effect of the SG massless effective mode, propagating at 
the speed of light, since dx/dt = 1 for the extreme points in Figure 5. Note also that the 
result (5.7) is consistent with our analysis in section 2.1, eq.(2.16). Let us study now the 
behaviour of the temperature corrections to this result in the DR, high-T limit. 

t-to 




Figure 3: Light cone regions for single charge 



t-t 




Figure 4: Light cone regions for two charges 



5.2 High-T corrections 

In section 4 we have discussed the DR high-T (and large distances) approximation to 
the retarded propagator. Let us write the contribution of the second term in the r.h.s of 
(4.18) to the current in (5.2) as 
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Figure 5: Current profile (first, free term), t > t + L 



A 2 2 d~(X) 



d/4 



(5.8) 



with given in (4.28). Only the cases h 2 and h 4 of (4.21) and (4.23) satisfy the physical 
constraint (4.20). Taking them in turn, h 2 of (4.21) gives 



~(X) 



3=1,2 



d 2 P c -ipo(t-to) jp(x-Xi) 1 



X 



p (po + «e) 2 - P 2 - Mo 



p 2 + m\ 



(Po + ief 



- 1 



(5.9) 



When the p integral in (5.8) is evaluated by the Residue Theorem for t > t , there 
is a contribution coming from the poles at p = ±oo p — ie and another one from the poles 
at p = —ie ± a/p 2 + ra\. The second one is always bounded. The remaining ^integrals 
can be expressed in terms of Bessel functions as 



A 



— E J2 SjSgnix-Xj) 



J=l,2 



9(t — \x — Xj\) j cLuJq 





rriT 



+ 9(\x — Xj\ — t) / duJo 
'o 



rriT 



Vt 2 - M 2 ] 

(5.10) 



for t = 0. The first term in the r.h.s of (5.10) renormalises the maximum value (5.5) 
for the free current. This is the effect of the renormalization of the residue of the free 
propagator discussed in section 4.3.3. As commented above, this gives rise to a current 
that in time-like regions rises to a constant independent of the length, i.e. a constant 
resistance. The second term in (5.10) is the effect of the "massive" mode my, yielding a 
causal oscillatory behaviour that is bounded in time. Our results are shown in Figure 6. 
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Figure 6: Time evolution of the second term — 5^ 1 'jb(X)/j max . The horizontal axis is mx 
and T = 5m, L = 1/m and A = 1. The curves in each plot are for times tm = 1.5, 5, 10 
and 20 with the later time plots extending further in mx. 

Finally, let us discuss the result of the calculation using (4.23) for the analytic con- 
tinuation of the propagator, which also satisfies the physical condition (4.27). Taking 
m s > rriT for definiteness, it is not difficult to see that the corresponding has three 
contributions. The first one is the pole at p\ = u 2 which gives again a renormalization 
of the factor multiplying 5^'j a which becomes 1 — \j^ r jm\ + {\ 2 / A^^/m 2 and therefore 
is of the same order as the first contribution in the r.h.s. of (5.10) (remember that m s 
and nif are of the same order and \3p -C m\) ). The second contribution comes from the 
poles at Po = p 2 + and gives exactly the same as the second and third terms in the 
r.h.s. of (5.10), multiplied by the factor (1 — m^/m 2 ) x ' 47r . Finally, the third contribution 
is given by the cut for p\ > p 2 + m 2 and gives 



71 rU S „•_■, n JO 



P 



5=1,2 

J^i (^; 2 -p 2 ) 2 (p 2 -B 2 + n^|./m^ , 1 ' 1 

for to = 0, where t' = m s t and x' = m s x. This latter contribution, because of the sin(A 2 /4) 
factor, is negligible in the small A limit compared to the other two. We have checked the 
equivalence between the results obtained with these two representations of the propagator, 
numerically for several values of the parameters in the DR limit. Therefore, we confirm 
explicitly that the difference between the physical transport coefficients calculated with 
different versions of the analytically continued propagator is negligible within the DR 
limit, as long as the condition (4.27) is enforced. 

The end result is that, once the transients have passed, 

6jM(X) = j max [0 n + 2(t-t )L- 1 e 23 + [(t-t )-(x-x 2 )}L- 1 6 21 

+ [(x-x 1 ) + (t-t )]L- 1 e 13 ] (5.12) 
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where 

with 6ki = 0\6f defined as before. For high enough temperature, <C m?, and so the 
effect of the non-leading terms in the temperature expansion is relatively small. This 
describes the regime in which chiral symmetry is approximately restored (the 'molecular 
phase' of [13, 10, 12]) and for that reason gives results close to those of the massless 
Thirring model. This condition is satisfied in all our plots. At fixed temperature the 
effect increases as the fermionic coupling becomes strong. 

In terms of the dual boson parameters the conductance G of the fermion field is now 



However, in terms of the fermionic mass m and coupling constant g 2 , as would follow 
from a diagrammatic description in terms of fermion loops, we find high temperature 
corrections to G of (2.17) of the form 



G ~ 1 U _ o (5 ^)/ ( ^) (i + ^ 2 A) 2 r 2 y^ 2)/( ^ 2) i 

~ 2tt(1 + g 2 /7r) I 1671- 2 \T 2 ) 



(5.15) 



The next nonleading term is down by powers of m/T. 

The power of duality in resumming series in g 2 that would be unobtainable otherwise 
is very striking, and is our main practical result. 



6 Conclusions 

We have analyzed the fermion conductivity at finite temperature in the massive Thirring/sine- 
Gordon (MT/SG) models. We have shown that the use of the dual degrees of freedom can 
avoid the usual problems associated with the infra-red behaviour of transport coefficients 
at finite temperature. In this case, the SG boson is the relevant quasi-particle mode in 
the limit of strong fermion coupling constant and the fermion conductivity can be related 
to the boson retarded propagator, which we have analyzed in detail. 

First, we have studied the imaginary-time SG boson propagator. Its leading order for 
small momentum at zero frequency has been related to the MT fermion charge density. 
Particular attention has been paid to the high T limit, for small boson coupling A (di- 
mensional reduction). In this limit, a thermal scale mr = A 2 T/2 emerges. This is the 
scale which screens the large distance behaviour of the propagator and therefore ensures 
that the results for the transport coefficients are well behaved in the IR limit. This is a 
common feature of high-T expansions, where the my scale comes from high-T diagrams 
of arbitrarily higher order in the coupling constant. In this sense, our results provide a 
partial resummation in A which ensure the IR finiteness. 

We have paid particular attention to the problem of the analytic continuation of the 
approximate (in DR) imaginary-time results. We have shown that by imposing the stan- 
dard mathematical conditions to the propagator only, one cannot guarantee that the 
difference between analytically continued propagators of approximately equal imaginary- 
time propagators remains small. However, one can ensure that physical results remain 
insensitive to this choice by imposing extra physical conditions. In particular, we have 
shown this in detail for the conductivity, where we demand that it remains bounded in 
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time. We have also discussed the role of the "thermal" modes in the retarded propaga- 
tor. Because of the ambiguities associated with the analytic continuation of approximate 
results, the position of poles and branch cuts in the retarded propagator cannot be de- 
termined with this approach. We have considered specific examples, where the analytic 
structure differs but the imaginary-time values are equivalent to leading order in DR. Our 
results are not in contradiction with considering uit a true thermal mass, but are also 
compatible with massless quasi-particles, since rriT is a soft quantity. The crucial point 
is that none of these interpretations prohibits us obtaining physically meaningful results 
for the conductivity. However, on the formal side, it is disappointing that we have had 
to apply to many approximations in what is an "exactly solvable" model. 

Once the physical condition is imposed, the fermion conductivity for this model can 
be studied perturbatively at high T. We have found that the resistance remains approxi- 
mately constant for long times, inside the causality region. The free SG field contribution 
correspond to an exactly constant resistance, while the high T corrections renormalize 
that constant and generate also bounded transient oscillations. 

Thus, we have managed to study in detail several aspects of this problem which 
will be generic to more realistic problems where the imaginary-time formalism is used. In 
particular, we have highlighted the importance of the analytic continuation of approximate 
results and the IR zero Matsubara mode. Studies of analytic continuation for propagators 
[30] or higher-order Green functions [34] are usually made in terms of exact functional 
forms and do not depend on the values at zero Matsubara energies. As we have noted, 
the IR sector is of vital importance to the physics and the only direct piece of information 
in this sector is in fact the zero Matsubara mode, the one not used in determining the 
analytic continuation. We have therefore proposed an additional condition to be used in 
these cases, (A. 30) which combined to physical conditions such as (4.27) allows us to find 
unambiguous physical predictions. 
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A Thermal propagators and Analytic Continuation 

The standard analysis is given in several places, such as [17]. We repeat it here to fix our 
notation and definitions, but also because we wish to go beyond the usual discussions to 
see how best to deal with approximate results. 

Take t e C, complex and lying on a given directed contour C with ends separated by 
—i/3. Define the thermal propagator as 

A c (x,t) = ((T c 0(x)0(O)» ee A>(x,t)9 c (t) + A<(x,t)6 c (-t), (A.l) 
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where A > (x, t) = (((j)(x) (f)(0))), A < (x,t) = ((0(O)0(x))), and 6c means time ordering 
along C as usual. The KMS conditions are then 

A > (x,t) = A < (x,t + if3). (A.2) 

The Fourier transforms of the propagators are defined to be 

A>(x,t) = j ^-J ge— * e ^A>(o;,p), (A.3) 

and similarly for other propagators. Then, KMS conditions in momentum space read 

A < (u,p) = e~^A > {u,p) (A.4) 

Now we introduce the spectral function 

p(cv,p) = A>,p) - A<(cv,p) = (1 - e"*") A>,p). (A.5) 

Note that in position space p(x,t) = ({[<f)(x), 0(0)])). On the other hand, A > (x, t) = 
A < (x, —t) so that p(x, —t) = —p(x,t) and then p(u>,p) = —p(—uj,p) so that 

=-p(o/,p)=0. (A.6) 

■oo 27r 

All versions of the propagator can be obtained from the spectral function. For in- 
stance, from (A.5) one readily has 

A>(u,p) = p(u,p)[l + f(u)}, 

A<(ou,p) = p(u,p)f(u), (A.7) 

where 

/(*) = ( A - 8 ) 

is the Bose-Einstein distribution. The simplest example is that of a free scalar field of 
mass m, for which 

Po(u,p) = 2-ne{u))5{u) 2 - p 2 -m 2 ). (A.9) 
More generally, let us introduce now the retarded and advanced propagators 

A R (x-x',t-t') = 9(t-lf){{[<!>(x,t),<j>(a/,{)])) (A.10) 
A A (x-x',t-t') = e(t-t){{[<f>(x , ,t),<i>(x,t)])) (A.ll) 

which we define for real time t—t' only because of the non-analytic nature of the Heaviside 
functions. 

Using the representation for the step function 

(A.12) 

with e — > + , the retarded propagator can also be written in terms of the spectral function 

as 

a , \ . f°° du' p(ou',p) 

A R (u;, P ) = i HK p > . A.13 

J_ 00 27T UJ - CO' + It 
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Now unlike the retarded function in real-time with its Heaviside functions, this energy- 
dependent retarded function has a simple extension into the complex energy plane. It is 
convenient to define a general function of complex energy 

This seems a little trivial as it simply related to both retarded and advanced functions 
through 

A R (u,p) = -iA(u + ie,p) 

A A (u,p) = iA(u-ie,p). (A. 15) 

with io e R. The point is that, on assuming that the grand canonical average ensures 
uniform convergence [34] of the thermal traces, we are assured that this generalised 
propagator function A(z,p) is bounded and analytic for all complex energies z except for 
real z, and A must tend to zero as \z\ — > oo [30, 34]. These are properties we will exploit 
below so A is a useful intermediate object to work with. 

Finally, let us introduce the imaginary time propagator Ay = A^ when C = [0, —i/3]. 
Then, changing variables to r = it G [0, j3\ C R gives 

A T (x,r) = A>(x, -it) = J ^-j ge— e^p(u,p)(l + f(co)). (A.16) 



Therefore, the Fourier transform of the IT propagator 

p/3 poo 

A T {iu n ,p)= dr dxe- ipx e luJnT A T (x,r) (A.17) 

J0 J-oo 

with uj n = 2nn/P, reads, from (A.16), 

f°° duJ> P(u',p) 
A T (iu n ,p) = - / — ^ ' PJ A.18 

It is crucial that we remember that in energy coordinates, the IT propagator is a sequence 
of functions of p and not actually a differentiable function of a continuous energy variable. 
This subtlety is easily missed since it is clear that we can relate the generalised propagator 
function, A, to the IT propagator through 

A(z = iu n ,p) = A T (iuj n ,p), neZ (A. 19) 

The key to analytic continuation is to note that we can not do the reverse easily, 
i.e. from the IT propagator we can not use this relationship alone to determine the 
full generalised propagator A at all complex energies. This is particularly confusing 
as analytic calculations, such as here, never give us a sequence of functions of p for a 
calculation of the IT propagator. Instead, we write down a function of a continuous 
complex variable, u n , and then note that for the IT propagator it is only to be taken 
at Matsubara energies. The obvious analytic continuation is to drop the restriction to 
Matsubara frequencies in such IT propagator expressions, and inspired by (A. 19) we 
would guess that we have then found the generalised propagator A. Unfortunately, this 
analytic continuation procedure is not unique as it stands. For instance, multiplying the 
IT propagator by arbitrary factors of exp(i/5u;„) = 1 would give different results for A 
but would not alter the values at Matsubara frequencies, so that (A. 19) holds. 

33 



There is, however, a way forward and we can find a scheme for analytically continuing 
from the IT function that gives the unique function A. The necessary and sufficient 
conditions were first stated by Baym and Mermin [30]. We simply quote here the result: 
given the IT discrete propagator A T (iu n ), then the unique function A(u) satisfies the 
following BM conditions: 

1. A(z = iuj n ) = A T (iou n ), \n\ e Z + 

2. lim {A(z,p)} = V lm(z) ^ 

\z\=oo 

3. A(z,p) is analytic V lm(z) ^ 0, tt (A.20) 

Therefore, if the result of an imaginary-time calculation can be written as Ax(z,p), 
analytic off the real axis and satisfying the second BM condition, it is then guaranteed 
to be the one generalised propagator function A(z,p). 

The retarded and advanced functions for real energy are then given by (A. 15), and 
a simple Fourier transform will then give the real-time retarded and advanced functions 
required for dynamical problems [17]. This is the usual situation when dealing with 
perturbative calculations, when the Matsubara sums in the loops are performed first [17]. 
However, the key point here is that we are dealing with nonperturbative expressions, like 
(4.4), which are nonperturbative in A, so that this procedure is not valid and finding the 
appropriate analytic continuation of the original imaginary-time expression is not an easy 
task. 

For instance, consider the free field case. The spectral function p of (A. 9), when 
inserted in (A. 18), gives the sequence of functions 

A° T (tuj n ,p)= 1 neZ (A.21) 

+ P + m 

One quickly sees that the function 

A°(z,p) = ^^ - (A.22) 

z z — p z — vn A 

obeys all the BM conditions (A.20) and therefore is the unique generalised propagator. 
The retarded propagator then follows from (A. 15) as 

A° R (u,p) = -——^- 2 -. (A.23) 

{uj + ley — p z — m A 

This is temperature independent, as it should since the dispersion law for free particles 
does not depend on the medium properties. 



A.l Analytic continuation of approximate results 

The results quoted above are all very well for a well defined function, such as the free 
propagator. However, in quantum field theory one does not have exact results for the 
Imaginary-Time Formalism (ITF) propagator at Matsubara frequencies. If we had been 
doing a numerical Monte Carlo calculation, the errors are partly statistical and random 
and partly come from the fact that only a finite number of Matsubara frequencies are 
calculated directly. This severely limits the accuracy of any possible analytic continuation 
from numerical data. 

In an analytic calculation the errors are usually of a functional form, that is we assume 
that the approximate result differs from the true result by an amount given by a suitable 
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analytic function. Indeed, we rarely calculate a sequence of numbers, but we invariably 
find a meromorphic functional form, even though we know it to be approximate and 
that it is strictly only valid for Matsubara frequencies. By checking that the function 
obeys the BM conditions (A. 20), we are then guaranteed that we have already found the 
relevant AC of our approximate IT propagator to the whole complex plane and hence the 
approximate retarded propagator. 

However, there is an interesting question one can ask. Suppose the functional error is 
order r\ where 77 may be a small coupling, a ratio of a mass scale to the high temperature, 
or whatever. Do the higher order, 0(rj), corrections at the Matsubara frequencies also 
lead to small corrections to its analytic continuation for all complex energies? In our case 
are the functional corrections at Matsubara frequencies, when analytically continued, still 
small at real energies? 

In fact the answer to this, posed in this simple minded way, is no. A small correction 
to the effective mass (real or complex) is a genuinely small change to the value of a 
propagator at all Matsubara frequencies. However, at real energies, near mass-shell, 
the propagator with and without the small shift will differ by large amounts. See our 
discussion in sections 4 and 5, or even the case of the free propagator analyzed below. 

Luckily though, this example shows that we are not interested in reproducing the 
analytically continued form of the function such that the errors in the value of the function 
are always bounded. Our goal is that the physics obtained from such functions be well 
represented. Our functions are Green functions that are not directly physical, but tools 
used to calculate scattering rates etc. It is not important if our Green function differs by 
an infinite amount from the true function at some real energies, provided the physical 
predictions are accurate. For example, our experience tells us that any error in our 
knowledge of the exact value of a mass leads only to comparatively small errors in cross 
sections. What mattered to the physics was that the propagator always had a single 
pole that is important and a small shift in its position is going to be a small error in 
physical results. In our case, the conductivity will turn out to pick up small 0(rj) errors 
even if these come from shifts in the mass value which leads to large differences in the 
propagators at certain values (near mass shell). 

We have highlighted this example as this is the task we face in our model. Namely 
we calculate an exact ITF propagator but to produce more understandable expressions 
need to approximate it, in our case by expanding in m/T. We need to be sure that 
terms dropped at Matsubara frequencies do not alter the physics which depends on the 
appropriate analytic continuation, not simply on the values at Matsubara frequencies. 

A. 1.1 The free case. 

To understand our approach, let us take the universal example of an exactly solvable 
quantum field theory, a free scalar field. We can compare the analytic continuation of 
approximations to the ITF version against the known exact result at real energies. Of 
course, in this case it seems a little stupid to make any approximation, such as DR, 
because we know the answer exactly. However, it will help us to see how we can be 
misled by DR. 

The ITF propagator is given in (A. 21). The unique generalised free propagator, 
the one obeying the BM conditions (A. 20), is just (A. 22). Suppose now that we were 
interested in the limit uj p <C 2-kT where u 2 = p 2 + m 2 (large distances and inverse masses 
compared to (3 or "soft, long wavelength" modes). This is what we will mean by high T 
here (for the full interacting theory we will also take T ^> m, the fermion mass, see main 
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text.). In that limit we find an approximate form for the propagator Ai, much as one 
has in actual calculations. This is given by 

A° T {iuj n ,p) « A%(iuj n ,p) 

i.e., this is the DR approximation and we pick up only the zero mode (n ^ represent 
"heavy" modes compared to uo p ). Now, if we find the unique generalised propagator, Ai, 
associated through (A. 20) to the IT function A^i, this analytic continuation of the DR 
propagator Ai is 

Ai(w,p)=0 forlmcj^0,7r =>- Af (u,p) = A? {uj,p) = 0, (A.25) 

since it is the behaviour of the IT function near the limit point at n — ±oo that controls 
the continuation, not the value at any one isolated finite n Matsubara energies. 

This is clearly not a very good approximation for the real-energy free propagator 
and is in fact inconsistent. First of all, we know that the true answer for the retarded 
propagator (A. 23) is T -independent. Thus, if we continue first and then take the high T 
limit, the answer is still (A. 23). Therefore, the DR and analytic continuation operations 
do not commute in this case. Second, even though (A. 24) is a good and well controlled 
approximation to the IT propagator at high T, that is not the case for its real-energy 
version in (A.25) when compared to the true answer in (A. 23). For any arbitrarily small 
uj, they differ by ou~ 2 , which is not small at high T. On the other hand for uj > 2nT the 
difference between the two is u~ 2 x O (u 2 /(2itT) 2 ) which is negligible to leading order 
in DR. Therefore, the conclusion is that by performing this DR style continuation, up to 
the accuracy of our approximations we are missing the right behaviour of the retarded 
propagator for frequencies uj < 2nT. 

Let us see if we can just solve the problem by making a better approximation. Consider 
A T2 where 

A%{iu n ,p) 
A T2 (iuj n ,p) 

This is an improvement over the simple DR form of (A. 24). Now all the modes appear 
as we take the leading high temperature term for each mode rather than working to 
a fixed accuracy for all terms. Notice however that the n = value does not fit the 
simple functional form we have given for the n ^ heavy modes. Thus we should not 
be surprised when the unique analytic function obeying the BM conditions (A. 20), when 
trying to match A T2 , does not fit the n = term. The continuation is controlled purely 
by the analytic behaviour of the Matsubara sequence near the limit points of ±ioo, and 
here we find 

A 2 (u,p) = -± (A.28) 
uj a 

The hard frequency behaviour is now a much better approximation to that of the exact 
answer, but it is still a bad approximation for the soft frequencies. This is not surprising 
as we have not used the one soft frequency value we know. 

The problem does not lie in the mathematical process of analytic continuation. Both 
(A.25) and (A.28) satisfy the BM conditions (A. 20) and they are therefore the unique 
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l + O 



ujz 



(2nTf 



, A° T1 (iuj n ,p) := ^ 



(A.24) 



= A T2 (iuj n ,p) 



l + O 



( "I 
\(2nT)< 



CO 2 ^ 



'nk 



(2nkTf 



(A.26) 
(A.27) 



continuations to the relevant generalised propagators of Ati and Ay 2 respectively. On 
the other hand the same is true for the exact answer (A. 23). Therefore, we have three 
apparently different continuations of the high temperature ITF propagator. This is not a 
contradiction, because we are dealing with approximate expressions for the IT propagator. 
Therefore, we are free to choose the continued version which gives a better approximation 
also for real frequencies. The problem is how can we specify the best continuation in a 
general case where we do not have the exact analytic answer? Is there any extra condition 
over and above those of Baym and Mermin (A. 20) which we can demand? 

The answer is that we should demand that the continuation also passes through the 
n = mode value, so incorporating the information about the n = mode. This value is 
unimportant for the continuations, since it is the behaviour of the ITF function near the 
limit points at ±ioo which controls them. However, the n = value is the only direct 
information we have on the soft energy region, so we ought to use it if we expect to get 
that physics correct. However, the continuation is fixed uniquely by the value of the 
functions at hard Matsubara frequencies and so we can not just choose its value at zero 
energy. What we must do is exploit the fact that we are working with approximations 
and alter the hard frequency values by small amounts equivalent to the inherent errors of 
the approximation. We will do this until we have a set of values at Matsubara frequencies 
which give an analytic continuation which agrees up to the errors of the approximation 
scheme at all Matsubara frequencies with the known values, including the zero energy 
one. 

Let's see how this works for the free propagator. For that purpose we recognise that 
adding a small energy independent constant will do the trick, but there is no unique 
choice. In this case one can see from the functional form for the ITF propagator what 
the best choice is. Thus, adding 'small' corrections in a high temperature regime to the 
form A2, we find a third form 



This is not a unique choice and here we have exploited the information in the ITF form. 
We therefore add a new condition for finding the unique generalised propagator through 
analytically continuing an IT propagator, namely 



Note that this condition demands that A R (0,p) = A A (0,p). However, this is guaranteed 
by the equal time commutation relations for two-point functions. An equivalent identity 
is true for higher-point functions [35]. 

In principle, all choices which differ by small amounts for each Matsubara value should 
give equivalent physical results, again up to the accuracy of the approximation scheme. 
However, as we have seen in section 4, in addition to the mathematical conditions we 
have just discussed, often one has to demand that the physical answer is meaningful. 
That is the case for the long-time behaviour of the conductivity discussed in section 5. 

B Algebraic analysis of the sine-Gordon propagator 

Here we will give more details about some of the calculations performed in section 3. 

First, we will show how to arrive to (3.12). Using (3.1) with the SG generating 
functional given in [15] and taking the /i — > + limit, we find 




(A.29) 



A(iu n ,p) « A- 



. T (iu n ,p) Vn>0. 



(A.30) 
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n=l 



(B.l) 



^2n{X\, . . . , X 2n ) 



A \ 2 / 1 x 2 



47T / V 71 ! 



2A 2 Vp 



x 



2n 



A 2 /4vr' 



2n 



2n 



2n 



n / ^ 



dou j i 



m=l 



^ln{[g 2 (^,r,)]^/ 47r } 

,z=i 
x 

(2n j-1 
j=2 fc=l 



^ln{[Q 2 (a; m -x,r m -r)r A / 4 n 



eje fc A 2 /47r 



(B.2) 



In order to obtain (3.12), we note that in the / and m sums in the above integral we 
can always relabel variables so that Xj becomes any of the x±, . . . ,x n if j < n and so 



on for Xj with n < j < 2n. The term Y\ 1 



<k<j<2n 



[Q 2 (Xj - X k ,Tj - 7fc)] 



€je k X 2 /4-n- 



remains 



invariant under such relabelling and also under the exchange Xj Xj +n , Vj = 1, . . . ,n. 
In this way, we end up with (3.12), where the InQ 2 factors are written only in terms of 
xi,x 2 and x 2n - 

Note that (B.2) is /xo independent. For convenience, in (3.12) we have written the final 
expression in terms of A again. They are equivalent in the no — > + limit. Let us focus 
now on the UV behaviour of (B.2), i.e, small x, r. This corresponds to the behaviour of 
the integrand in near the regions where Q 2 — > + in denominators and logs (see a similar 
discussion in [15] for the partition function). On the other hand, Q 2 (x, r) ss 7r 2 T 2 (:r 2 + T 2 ) 
for (x, t) — > (0,0). Thus, clearly the most divergent contributions arise when x = r = 
and come from integrals like: 



/ d2y f 

JT,y~0 JT,y'r 



2„/ 



dfy 

v [Q 2 (y-y')} 



A 2 /4tt 



\nQ 2 (y)\nQ 2 (y>) 



In 2 r 



(B.3) 



where J T d 2 y = J Q cfa/ f^dyi and r 2 = y 2 + y 2 . Clearly, the above integral converges 
near r = for A 2 < An. 

Next, we will give more details on the derivation of (3.25). When taking the Fourier 
transform of (B.2) we realise that the only dependence with (r, x) appears in ^fc=i e fe m Q 2 { x k~ 
x, Tk — t) = — AnY^k=i e k^r( x k — x , T k — t), so that we can write (3.21) with 



A 2 Z»(T) 



^ \n\ 



a ( T 



2A 2 \p 



A 2 /4tt' 



2n 



^(^n,p), (B.4) 



and 



2" pP r+oo 

E^(iuj m ,p)A^(iuj m ,p) = IJ/ dr J / dx 3 HI [Q 2 ^ - ^fc^i - r fc)] 

j=l "'-° l<fc<j<2n 

2n 



eje fc A 2 /47r 



X 



J2e j \nQ 2 (x j ,T j ) 



2ra 



k=i 



(B.5) 



Let us now write the expression (B.5) in a slightly different form, more convenient for 
our purposes. First, we relabel the variables in the j and k sums as we have just done 
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with (B.2). Thus, separating the terms with j = k and j ^ k in the sum, we can write: 

2n „fs 



/p i>+oo ^ 

oItj dXj Yl [Q 2 ( X j ~Xk,Tj - Tk)] 

l<k<j<2n 



InQ (xi,ri) 



x {ne^^e"^ 1 + n(n - l) e ^ mr2 e - ip:C2 - n 2 'e iu}mT2n e ~ ipX2n .} 



Now let us change variables to: 



so that 



(B.6) 



(B.7) 



E x 'i ' x i ~ Xk = E ^ ( k< j< 2n )' 

i=fc+i 



and so on for r^, . . . t[. Note that with this change of variable, 



(B.8) 



l<k<j<2n l<k<j<2n 



Q 2 ( E ^ E ^ 

;=fc+i l=k+l 



e j e fc A 2 /4vr 



(B.9) 



is independent of (x[,t[). This is just a consequence of two-dimensional translation 
invariance or, equivalently, total energy-momentum conservation in any diagram, which 
becomes manifest with this change of variables. Thus, we can separate the free propagator 
A^(iu m ,p) in the r.h.s of (B.6) yielding (3.25). 
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